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Cavitation in holographic sQGP 
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We study the possibility of cavitation in the nonconformal M — 2* SU(N) theory which is a 
mass deformation of J\f = 4 SU (N) Yang- Mills theory. The second order transport coefficients are 
known from the numerical work using AdS/CFT by Buchel and collaborators. Using these and the 
approach of Rajagopal and Tripuraneni, we investigate the flow equations in a 1+1 dimensional 
boost invariant set up. We find that the string theory model does not exhibit cavitation before 
phase transition is reached. We give a semi-analytic explanation of this finding. 



I. INTRODUCTION 

Viscous relativistic hydrodynamics has been applied 
successfully to the strongly coupled Quark Gluon Plasma 
(sQGP) created in nuclear collisions at Relativistic Heavy 
Ions Collider (RHIC). The fact that a low shear viscosity 
to entropy density ratio seems to be favoured from this 
experiment coupled with the fact that any holographic 
transport calculation using AdS/CFT also gives rise to a 
low value for this ratio [1, 2] has triggered a huge amount 
of activity in using AdS/CFT methods to understand 
cousins of the sQGP (see [3] for a recent review). Since 
then, other hydrodynamic coefficients have been obtained 
using AdS/CFT. Of particular interest for this work is 
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the bulk viscosity £ and the relaxation times t% , Tj 

While the bulk viscosity is usually very small in sQGP, 
there is good evidence that it rises sharply near the criti- 
cal temperature [4-8] (see also [9, 10]). With high enough 
viscosity, the boost-invariant solution to the hydrody- 
namics equations become unstable against small pertur- 
bations [11]. A related effect arises due to the fact that 
the bulk viscosity affects the trace of the stress-energy 
tensor and it can lower the effective longitudinal pressure. 
The longitudinal pressure of the plasma could become 
negative making the plasma unstable to clustering and 
cavitation. The proposition is then that this cavitation 
process could be one possible description of hadroniza- 
tion of the sQGP. 

Cavitation occurs when a fluid has such a large gradi- 
ent velocity and bulk/shear viscosity coefficient that the 
viscous correction to the pressure is large and the total 
effective pressure dips below the vapour pressure of the 
fluid. At this point, bubbles of the gaseous phase are 
formed and the hydrodynamics description fails. Cavita- 
tion is a well studied phenomena that occurs, for exam- 
ple, in pumps and propellers where the induced shock- 
waves are usually harmful. It has numerous technological 
application ranging from shock wave lithotripsy (a kid- 
ney stone treatment) to ultrasonic cleaning baths. 

While cavitation is fairly well understood, it is not de- 
scribed by hydrodynamics which is based on an expan- 
sion in gradient of velocity. When the viscous terms are 
of the same size as the perfect fluid terms we presumably 
lose control of the perturbative expansion in gradients. 



In the context of RHIC, as the plasma expands it cools 
until a crossover occurs at a critical temperature which 
we take to be T c — 190 MeV. Finding that a fluid cav- 
itates is therefore the same as figuring out whether the 
hydrodynamics description breaks down before we reach 
the critical temperature 1 . 

This question has been looked into (with RHIC intial 
data in mind) using a boost invariant longitudinal flow 
(1 + 1 dimensional) in [13-15] (see [16] for an analysis in 
the context of neutron stars) . The analysis was extended 
to include transverse flow in a 2 + 1 dimensional code by 
Song and Heinz [19]. In this paper, we will follow the 
(1+1) dimensional analysis of Rajagopal and Tripuraneni 
(RT) [15]. Their approach uses a mixture of inputs from 
lattice simulation of QCD as well as results borrowed 
from the conformal plasma: 

1. They use a phenomenological model for the equa- 
tion of state fitted to lattice data and they assume 
a vanishing chemical potential. 

2. They assume a phenomenological bulk viscosity 
(again from lattice calculations) which rises sharply 
near the crossover temperature while the shear vis- 
cosity to entropy ratio is assumed to be a few fac- 
tors within the leading order AdS/CFT result of 

1/4-7T. 

3. For the second order transport coefficients (like re- 
laxation times) they take the AdS/CFT values for 
a conformed plasma. 

They find that the rise in bulk viscosity can lead to 
cavitation while pointing out that large relaxation times 
could hinder this effect. The analysis of RT will be the 
benchmark to which we will compare in this paper and 
we denote their results QCD — RT in the various figures. 
It would be interesting to study cavitation in a theory 
close to QCD where the AdS/CFT let us calculate and 
predict all the relevant coefficients. 



1 Assuming cavitation does occur, it could have important phe- 
nomenological implications on photon and dilepton production 
[12]. 
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In this paper, we look at cavitation in Af = 2* model 
whose transport properties have been investigated in 
some details by Buchel et al. [17] and Buchel and Pag- 
nutti [18]. In particular, we study the 1+1 d boost invari- 
ant expansion following Rajagopal and Tripuraneni but 
with the transport parameters as dictated by this model 2 . 
This model mimics several properties that the sQGP is 
thought to have. It shows a rising bulk viscosity near 
the phase transition temperature and a low shear viscos- 
ity to entropy density ratio. While these features are in 
common with [15], the relaxation times are not-they rise 
sharply near the phase transition temperature, quite un- 
like what happens in the conformal case where they are 
inversely proportional to temperature. This behaviour, 
which has been anticipated by Song and Heinz [19], is re- 
sponsible for a "critical slow down" . We find that using 
a normalized data for Af = 2* which sets T c « 190 MeV 
and using an initialization time of Tj ~ 0.5/m/c, cavita- 
tion does not occur in this stringy example. The main 
reason for this is the lower value of bulk viscosity and 
higher value of pressure in these models as compared to 
the ones used in RT. 

We then repeat the analysis of RT using a phenomeno- 
logical formula for the relaxation time that exhibit criti- 
cal slow-down as seen in the N — 2* plasma 
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T cV /l - T c /T 



(1) 



The Af = 2* relaxation is well approximated by the for- 
mula above (near T c ) with C — 0.08. We find that cavi- 
tation is prevented for values of C > 2. We also analyze 
the situation where the equation of state is taken from 
the lattice but the bulk viscosity and the relaxation time 
are taken from the Af = 2* data. We find that cavitation 
is prevented mainly because of the lower bulk viscosity 
of the Af = 2* compared to the lattice result. 

Our conclusion is that critical slow-down hinders cav- 
itation but one needs to have a bigger relaxation time 
than the naive value given by the Af = 2* toy model. It 
is interesting that there is no range of parameters that 
makes the toy model Af = 2* cavitate before it reaches 
the phase transition. We give a semi-analytical proof of 
this behavior. 

This paper is organized as follows. We review the 
second order hydrodynamics equations in §11. In §111, 
we consider the possibility of cavitation in a conformal 
plasma and give an analytic argument that this cannot 
happen. In §IV, we turn to the numerical study of the 
1+1 d boost invariant flow equations using the Af = 2* 
model. We conclude with a discussion in §V. The basic 
assumptions used in RT are reviewed in an appendix. 



2 The numerical data has been kindly supplied to us by Alex 
Buchel. 



II. BOOST INVARIANT VISCOUS FLUID 

We are interested in a relativistic fluid whose stress- 
energy momentum tensor is 



(2) 



where e and p are the fluid energy density and pressure 
respectively, w M is the fluid 4-velocity (with u^u^ = 1), 
IP" is the viscosity tensor and A^ v = g^ v - u^u v is the 
projection tensor orthogonal to the 4-velocity. The vis- 
cosity tensor is also orthogonal to u M with u^IP" = and 
g^ v is the metric with signature (+,—,—,—). The equa- 
tion of motion of hydrodynamics is just the conservation 
of stress-energy 







(3) 



where V M is the geometric covariant derivative. We in- 
troduce the following notation 

D = u^, , (4) 

A <ttv> = ^(A^A 1 ^ + A^ 3 A va ) - ^A^A^ A a p . 

We will be considering a purely longitudinal flow. The 
conservation of stress-energy in the longitudinal direction 



Ufl T^ = is 



De + (e + p) V • u + n""V M u„ = 



(5) 



The most general viscous tensor at second order in per- 
turbation theory has been worked out in [20] . Working in 
flat space with no vorticity reduces the number of terms. 
Writing IP" = ir^ v + A^"TI we find 

= 7 7 V < ^^ > - tjt^ ^D(V<V>) + ^pV^u 1 ^ 
-X 1 V <t "u x V > u x (6) 



while the trace is given by 

n = -cv • u + Ct£d(v • u) 



(7) 



where in the last expression we neglected two second or- 
der corrections which are not zero even in our purely 
longitudinal flow (the coefficient £i, 2 in [20] are set to 0). 
This is a common procedure. The coefficients in front of 
each viscous terms must be determined experimentally or 
via some theoretical model. They can depend on temper- 
ature or time and they are thus general scalar functions. 
They are the shear viscosity 77, the bulk viscosity £, the 
relaxation times t^,t^ and Ai. 

At this point this second order hydrodynamic (just like 
the first order) is acausal. The viscous tensor responds 
instantaneously to change in the local fluid properties. It 
has been shown that the superluminal modes are not fun- 
damentally a problem since they corresponds to modes 
outside of the regime of the effective theory of hydrody- 
namics [21-23]. Nevertheless they pose problems when 
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one tries to solve the system of equation numerically. 
Mueller, Israel and Stewart [24-27] have proposed a so- 
lution which consist of extending the set of variables to 
include II and n 1 *" and transforming their constitutive 
relations Eq. (6) and (7) into differential equations. 

To do this, we adopt the approach of [28] (see also 
[29]) where we simply put the first order solution in the 
second order terms. The first order solutions are Tt^ u — 
r]V <fJ, u" > and II = (\7 • u and substituting in the second 
order terms, 

2 A ' 



n 



-CV • u - t^DTI 



(8) 



where we have used the relation Drj = —77V • u. This 
follows from assuming that rj oc s where s is the entropy 
density which obeys Dins = —V • u [20, 28]. In the 
second equation, we have dropped the term — t^V • uDC, 
which we checked had a very small effect on our analysis. 
Eq. (5) and (8) now form an hyperbolic set of equations 
which is causal and well posed for numerical methods. 

As we already mentioned, we are interested in a boost 
invariant longitudinal flow in the z direction. Instead of 
Minkowski coordinate (i, z), it is convenient to change 
to (t, £) where r = \Jt 2 — z 2 is the proper time and 
£ = arctanh(z/t). In these coordinates the metric is 
g^v = diag(l, -1, -1, -t 2 ), = (1,0,0,0) and boost 
invariance implies that all quantities are independent of 
£. Assuming boost invariance and longitudinal flow sim- 
plifies the stress-energy tensor to 
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(n-*)/rV 



It will be useful to define the transverse and longitudinal 
effective pressure as 

1 



p± = P + n + ^$ , (10) 
p e = p + n-$ . (11) 

In this coordinate system, the evolution equations are 

p + Il - $ 



de 

dr t 
n dr 3r 3t 2t7 2 



dr 



r 



(12) 
(13) 
(14) 



III. CAVITATION IN A CONFORMAL PLASMA 

To get an idea of what we are going after, let us start 
with a conformal plasma (£ = and p — e/3). Consider 
the situation where the hypothetical cavitation occurs at 
time t = Tf corresponding to temperature T — Tf and 
substitute 

/Q rr.4 V b c ds 

in Eq. (12) and (13) where b,c,d are unspecified numer- 
ical coefficient and s is the entropy density. Consider a 
series solution for T and $ around r = Tf. This is given 

by 

T = Tf -?L(T-Tf) + ~ • , $ = e jTf + cP(T-T f ) + - ■ ■ , 

(15) 

with 



9b 2 C1T 2 Tf 



The longitudinal pressure p — $ 

eoT 4 



[&V(4c + 3T f T f ) + 6dir 4 T f T f - 46 j tt)] . 

(16) 



p — $ = (r — Tf) 



9b 2 C1T 2 Tf 



[-46 3 tt 



+b 2 TT 2 (c + 3T f T f ) + 6dn i T f Tf] , (17) 

so that in order to get cavitation we need for t > Tf, 
p — $ < 0. This translates into 

6 2 7r(c7r - 46 + 3T f Tfir) + Qd^TjTf < . (18) 

From this we can conclude the following. 

• Increasing rj/s aids in cavitation. 

• Increasing t% hinders cavitation. 

• Increasing Ai hinders cavitation. 

Now the last term proportional to d is positive so it is a 
necessary condition that the term proportional to b must 
be negative for cavitation to occur. There will not be 
any cavitation unless 



1 3tt 

> 



TfTf 46 — 7TC 



(19) 



Now imagine that after r = r, the system has a hydro- 
dynamic evolution until cavitation occur at r = Tf. The 
cavitation temperature should be greater than the phase 
transition (crossover in the QCD plasma case) tempera- 
ture Tf > T c , 



1 1 

r = — — < 



(20) 



Now at RHIC, T c - 200 MeV, n - 0.5/m/c = 
1 / 400 MeV such that the bound roughly translates into 



r < 2. 



(21) 
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Plugging in the conformal values 6=1 and c = 2 ^ 1 ^ 2 ^ , 
we have r ~ 2.8 so we conclude that this condition can- 
not be met and cavitation does not occur. In fact since 
the term proportional to d is positive and comparable to 
the rest of the terms, its presence only makes cavitation 
harder. 

We now turn to the analysis in QCD and Af = 2* where 
we know the bulk viscosity is not zero and cavitation is 
more likely to happen. 



QCD-RT 



IV. NON-CONFORMAL CASE 

We compare the behavior of Af = 2* to the results of 
RT (which we denote QCD — RT in the graphs). Let us 
summarize what they did (see appendix A for details). 
They took a phenomenological approach where all pa- 
rameters have their conformal values with the exception 
of the bulk viscosity and the equation of state. They then 
scan the parameter space. Clearly the presence of bulk 
viscosity enhances the chances of cavitation since at first 
order II = — C/ r an d Pt oc II. The relaxation time on the 
other hand tends to slow the rise of IT and it therefore 
hinders cavitation. 

They find that cavitation indeed sets in before T c . The 
weakest link in this analysis is the choice of the relaxation 
time which was imported from the conformal plasma 
calculated using AdS/CFT. There exists a well known 
AdS/CFT model studied extensively by Buchel and col- 
laborators called the Af = 2* model which has a bulk 
viscosity [17, 18]. It is interesting to ask whether this 
model exhibits any sign of cavitation along the lines of 
Rajagopal and Tripuraneni's analysis. Let us first review 
some salient features of this model. 



Af = 2* 
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FIG. 1: Comparison between the M — 2* bulk viscosity and 
the phenomenological bulk viscosity fitted to lattice QCD cal- 
culation and used in RT 
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FIG. 2: Comparison between effective relaxation times. The 
Af = 2* model shows a sharp rise near T — T c , a sign of 
critical slow-down. 



Af = 2* is a mass deformation of = 4 and is dual to 
the Pilch- Warner geometry in IIB. Typically one needs 
to resort to numerics as the resulting equations are too 
hard to solve exactly. On the gravity side, there are 
two scalar fields in a certain potential. One can give 
masses mf,,m/ to bosonic and fcrmionic components of 
the Af = 2 hypermultiplet respectively. We will consider 
nrif — with TTi 7^ 0. When m 2 < m| this model ex- 
hibits a second order phase transition, in our case when 
mb/T c — 2.33. We will normalize the free energy such 
that the phase transition occurs at 190 MeV and choose 
the number of colour to be JV C = 3. At finite temper- 
ature, the fluid/gravity correspondence analysis of this 
model yields a bulk viscosity which is thought to respect 
the Buchel bound (/rj > 2(1/3 — c 2 s ) where c s is the ve- 
locity of sound [30]. The bulk viscosity rises up sharply 
near the phase transition (but not as sharply as in QCD). 
This is illustrated in Fig. 1. 

Another crucial difference is the behaviour of the re- 
laxation times as a function of T/T c . The Af — 2* model 



predicts the following effective relaxation time 



T eff 



3C r C 



1 



3 c 

4 n 



which exhibits a peak 

r^fT^ll-TjT]- 1 / 2 



(22) 



(23) 



near T — T c . Note that from the gravity side of the Af = 
2*, we can only obtain this particular combination of the 
two relaxation times. We parametrize our ignorance of 
the repartition between the two times with the parameter 

7, 



1 



1+ 7 ii 



3C 
4 rj 

H Teff 



r n =7 



1 



1 + 'T— — 



3( 

In 

IZ Teff 



(24) 







For 7 = 1, we have that = r^, for 7 = 00, t% 
and t£ = r e ff + lj while for 7 = 0, we have r£ = 



5 



and t% = r e ff fl + In the following we will test 

these three different values of 7. The behaviour of r e ff 
as a function of temperature is illustrated in Fig. 2. As 
we see, the effective relaxation time rises sharply as T c 
is approached. That this plays an important role was 
already anticipated in the work of Song and Heinz [19]. 
Together with the fact that the bulk viscosity in this 
model does not become too big near T c with (/r) 0.8 
makes cavitation harder. 

As in RT we choose the following initial conditions 
3>(t,-) = 0,n(ri) — 0, Ti — 0.5 fm/c and we take the 
initial energy density to be e = 15.8 GeV/fm 3 . This 
initial energy is a bit uncertain but we find that it does 
not affect the answer greatly. For the lattice equation of 
state Eq. (Al), this initial energy density corresponds to 
an initial temperature of T = 305 MeV. For the Af = 
2* equation of state, the initial temperature is lower at 
T = 245 MeV. Given that we stop the evolution at the 
same temperature in both cases (when T = 190 MeV), 
the evolution of the Af — 2* is shorter. 

In Fig. 3, we plot the longitudinal pressure for the Af = 
2* plasma compared to RT. The longitudinal pressure 
decreases with time but remains positive all the way up to 
T c . We have checked that the reason for this is primarily 
the smaller values of bulk viscosity in this model. In 
fact substituting the Af = 2* relaxation times into the 
analysis of RT still leads to cavitation in their model. In 
Fig. 4, we show how the longitudinal pressure is affected 
by varying 7 between the three limiting values of 0, 1, 00 
for the Af = 2* model. For 7 = 00, the behavior is quite 
different but the system still does not cavitate. 
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FIG. 3: The longitudinal pressure for Af = 2* compared with 
the findings of RT. The string model reaches the critical tem- 
perature faster and does not cavitate before that. For the 
QCD-RT model, cavitation occurs near r ~ 2.5 fm/c 




FIG. 4: The longitudinal pressure for N = 2* for different 
values of the parameter 7. 



some reasonable approximations, we will be able to de- 
rive a simple inequality for the non-conformal Af = 2* 
plasma for cavitation to occur. Let us assume that cavi- 
tation occurs at time r = Ty such that T > T c . We will 
further assume that the plasma is close to conformality 
or in other words p ~ e/3 which is accurate within 10% 
error for the Af = 2* plasma but is substantially different 
for the QCD plasma. For the QCD plasma considered in 
RT, the pressure works out to be smaller compared to the 
energy and it turns out that it will make the inequality 
only stronger (in other words cavitation occurs for the 
QCD plasma of [15] at an earlier time than what will be 
predicted by this analysis). We begin by substituting 



77 b 

S 47T 



P = e/3, e = e T 4 
We take Ai = and 7=1 for simplicity and 



n C 



r - :=r no + T ni( T - r /)' C = Co + Ci(r - Tf) , 



in Eq. (12) , (13) and (14). Consider a series solution for 
T , <f> and II around r = Tt. This is given by 



T 



2(4^) 
l2e Tf 



and 



n = (ci 



eo)-^ +7ri ( r ~ r /) + 



-*-/)+• 
(25) 

(26) 



Here we have parametrized the deviation from confor- 
mality by ei — €q. Now we solve for <fi and tti and we find 
that the longitudinal pressure around £q ~ e i is 



B. No cavitation for Af = 2* 

As in the conformal case where a series solution en- 
abled us to derive an inequality on the parameters, with 



Pt = 
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ir(3T Tfe 
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( T - T /)- 
(27) 
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Now the analysis is exactly as in the conformal case. In 
this case defining (with e/ = c\t=t { ) 



1 



T 



> 



3ne 



f 



f T f 



4be f + % nT f - e/Trr^I) 



(28) 



and choosing the RHIC initial time and temperature as 
in the conformal case leads to cavitation happening only 
if T£ < 2. Some comments are in order regarding the 
choice of ei k to. Our numerical analysis show that the 
TV = 2* plasma has II — > at late time corresponding 
to ei = €q. This is both because of lower bulk viscosity 
and because of large relaxation time. In RT, on the other 
hand we have IT ss — $ at late time (so ei ~ £o/2). So 
our current approximation is only good for the Af = 2* 
plasma. As we see in Fig. 5, cavitation does not happen 
for Af = 2* while RT on the other hand clearly cavitates. 
Since the QCD plasma of [15] is far from conformality, 
our approximation is rather crude and predicts cavita- 
tion to occur when T is between 1.3T C and 1.4T C while 
the numerics gives around 1.7T C . The reasons for this 
discrepancy is that the QCD pressure is lower than the 
conformal approximation and t\ ~ eo at late time. Ac- 
counting for both effects lead to an earlier cavitation time 
for RT. 
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FIG. 5: as a function of temperature. Cavitation occurs 
when r( < 2. 



C. Critical slow-down 

It is interesting to ask if the critical slow down phe- 
nomena where the relaxation time shoots up near T c can 
inhibit cavitation in the QCD plasma considered in RT. 
We take the equation of state and bulk viscosity from 
the lattice simulation as in Eq. (Al) and (A3) but in- 
stead consider a phenomenological relaxation time of the 
form [17] 



C 



Wl - T c /T 



(29) 



For C = 0.08 we find that this is a good fit to the Af = 
2* data for a large range of temperature although it is 
important to mention that the relaxation time does not 
actually diverge in the Af — 2* plasma at T — T c . As 
can be seen from Fig. 6, cavitation is prevented when 
one reaches value of C > 2. We conclude that critical 
slow-down does prevent cavitation but one needs a bigger 
effect than what is seen in the Af — 2* plasma. 

GeV 



P f [ 1 




FIG. 6: This is the longitudinal pressure for a plasma with 
equation of state Eq. (Al) and (A3) but with a phenomeno- 
logical relaxation time that exhibit critical slow-down for dif- 
ferent value of C. 



V. CONCLUSION 

In this paper we investigated if the Af — 2* model in 
string theory exhibits cavitation and we contrasted our 
finding with the analysis in QCD following closely [15] 
(RT). We found that for this choice of normalizations 
and initialization conditions this model docs not cavitate 
before T c is reached. Primarily, the smaller values of bulk 
viscosity and secondarily, the sharp rise in the relaxation 
time near T c are responsible for this behaviour. It is also 
likely that quantum effects as in [31] which correct rj/s 
and £ will play an important role in these studies. 

If we made the initial time Ti of hydrodynamic evo- 
lution smaller, the chance of cavitation becomes higher. 
For example, for the Af — 4 plasma, our analysis shows 
that cavitation could occur for an initialization time less 
then 0.3 fm/c while the Af = 2* plasma could cavitate for 
Ti < 0.14 fm/c. Simulations at RHIC suggest Ti closer to 
1 fm/c in which case cavitation would not happen. 

As stated in the introduction in order for the longitu- 
dinal pressure to become negative, a second order con- 
tribution to the pressure has to be as big as a zeroth 
order contribution. As such, mathematically the gra- 
dient expansion approximation used becomes incorrect. 
The final value of in our simulations (see Fig. 3) was 
around 20% of the initial value indicating that the nu- 
merical approximations starts to become suspect. This 
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raises the question if the < condition is a good 
marker for the breakdown of the hydrodynamics equa- 
tions. One can perhaps answer this question more accu- 
rately using AdS/CFT where a breakdown of the hydro- 
dynamics equations is signalled by the existence of some 
instability in the gravitational theory. 
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Appendix A: QCD equation of state 

Here we summarize the important assumptions in RT 
[15] that we used as a baseline for comparison throughout 
our analysis. The equation of state is taken from the lat- 
tice calculation of [32] and it can be nicely parametrized 



as follows 



3p 



rp4 



1 - 



l (l+cx P (^))' 



with d 2 = 0.24 GeV 2 , d A = 0.0054 GcV 4 , Ci = 
0.2073 GeV and c 2 = 0.0172 GeV. The crossover tem- 
perature is chosen to be T c = 190 MeV. We also have 



P(T) 

rp4 



dT'- 



3p 



To 



y/5 



(A2) 



For the various transport coefficients, RT choose 
„ 2-log2 



1 

47T 



Ai = 



2irT ' '"■ 8?r 2 T 
The bulk viscosity is taken from the lattice [33] 

2 



- = aicxp [—^7- 



0,2 



(A3) 



with ai = 0.901, a 2 = 0.061 and AT = RT varied 

their parameters and found the following: 

• Cavitation is insensitive to changes in t% and Ai 

• Cavitation disappears when is increased by 40 
from its baseline value taken to be = t%. 

The latter is in agreement with our finding that cavi- 
tation does not occur if we allow for a critical slow-down 
of the form Eq. 1 with C > 2. 
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